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This paper proves that non-convex quadratically constrained quadratic programs can be solved 
in polynomial time when their underlying graph is acyclic, provided the constraints satisfy a certain 
technical condition. When this condition is not satisfied, we propose a heuristic to obtain a feasible 
point. We demonstrate this approach on optimal power flow problems over radial networks. 
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Index Terms 

Optimal power flow, distribution circuits, radial networks, energy storage, SDP relaxation, minimum 
semidefinite rank. 

I. Introduction 



A quadratically constrained quadratic program (QCQP) is an optimization problem in which 
the objective function and the constraints are quadratic. In general, QCQPs are non-convex, and 
therefore lack computationally efficient solution methods. Many engineering problems including 
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optimal power flow (OPF) can be represented as QCQPs with complex variables. The contribution 
of this paper is to expand the class of non-convex QCQPs for which globally optimal solutions 
can be guaranteed. 

There is a large literature on optimal or approximate algorithms for QCQPs. One such method 
employs a convex semidefinite program that is a rank relaxation of the given QCQP This is 
commonly referred to as semidefinite relaxation (SDR). Such semidefinite programs are solvable 
in polynomial time using interior-point methods |TJ-[|3j. In some instances, an optimal solution 
of the original QCQP can be recovered from an optimal solution of its SDR. In other cases, SDR 
provides a way to approximate the solution of a QCQP. Thus, SDR provides a computationally 
tractable way to approach QCQPs Q, [|5j. For example, SDR has been applied to a variety of 
engineering problems such as MIMO antenna beam-forming J6J-J9J, sensor network localization 



[10], principal component analysis JTTJ and stability analysis p2[ . SDR has also been extensively 
used in systems and control theory applications p"3| , p"4| . Several authors have investigated exact 
relaxations, e.g., p"5|, flT6|, while others have applied SDR-based approximation techniques to 



NP-hard combinatorial problems and non-convex QCQPs, e.g., [ 17|-[ 19 1. The accuracy of these 



approximations has also been extensively studied, e.g., [20|-[22|. 

In this paper, we prove a sufficient condition under which QCQPs with underlying acyclic 
graph structures admit an efficient polynomial time solution using its SDR. We then apply 
our result to the optimal power flow (OPF) problem on radial networks. OPF is generally 
a non-convex, NP-hard problem that seeks to minimize some cost function, such as power 
loss, generation cost and/or user utilities, subject to engineering constraints. Since the original 
formulation of Carpentier in 1962 [23], various solution techniques have been used for this 



problem; we refer the reader to [24|-[28| for some surveys. OPF can be cast as a QCQP. The 
authors in [29|, [ 30 1 propose to solve its SDR as an approximation. The authors in J3TJ, instead, 
propose to solve its (convex) Lagrangian dual and provide a sufficient condition under which 
an optimal solution of the OPF can be recovered from a dual optimal solution. For IEEE test 
systems and other randomly generated circuits, these approaches have been shown to solve OPF 
optimally. Recently, OPF over radial networks has been of considerable interest. A checkable 



sufficient condition has been proved in [32|-[34| where an optimal solution of OPF can be 



recovered from an optimal solution of its SDR. This paper extends the previously known class 
of OPF problems that can be solved efficiently. 
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The paper is organized as follows. In Section [TTJ, we prove a sufficient condition for a non- 



convex QCQP over acyclic graphs to be solvable in polynomial time. In Section III we apply this 



result to characterize a class of OPF problems over radial networks that can be solved efficiently. 



In Section [TV] we describe a heuristic method to obtain feasible solutions for QCQPs that do not 
meet these conditions and thus their optimum solution cannot be directly recovered by solving its 
SDR. We further apply this technique for the OPF problem and demonstrate through simulations 
that we can always find a near-optimal feasible point for OPF. We conclude in Section |VJ 

II. QCQP'S AND SEMIDEFINITE RELAXATION 

Consider the following QCQP with complex variable x E C n . 
Primal problem P: 

minimize x H Cx 

xec n 

subject to: x n CkX < bk, k E K. 

where x n denotes the conjugate transpose of x, C is either an n x n complex positive definite 
matrix (denoted as C y 0) or a positive semidefmite matrix (denoted as C >z 0), /C is a finite 
index set, Ck is an n x n complex Hermitian matrix and bk is a scalar for each k E K.. 

If the matrices Ck, k E fC are positive semidefmite, then problem P is a convex program 



and can be solved in polynomial time [35|, [36 1. If, however, these matrices are not necessarily 



positive semidefmite, P is non-convex and NP-hard in general. The main result of this paper is 
the identification of a class of QCQPs that can be solved in polynomial time even though the 
matrices Ck, k E K, are not necessarily positive semidefmite. This result is applied in the next 
section to the optimal power flow problem on radial electric networks. 

We begin with some notation. Let [n] := {1, 2, ... , n} and for any matrix H, let Hy represent 
the element corresponding to the i th row and the j th column. Define a function Q from QCQP 
problems to undirected graphs as follows. For a QCQP problem P, the undirected graph Q(P) 
has vertex set [n] and the edge set defined as follows: 

(i, j) is an edge in Q(P) ^=^ (i ^ j) and (CV,- ^ or [Ck]ij ^ for some k E K) (1) 

Since the matrices C and Ck are Hermitian, the edges of the graph are undirected. Also, Q(P) 
has no self-loops from a vertex to itself. We restrict attention to QCQP problems P for which 
the graph G(P) is a tree, i.e., it is connected and acyclic. 
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For any vector of real numbers a, let a ^> denote that all elements of a are strictly positive. 
For any set of complex numbers U = {ui, . . . ,u r }, the relative interior of the convex hull J36| 
of the set is defined as: 



r 



{a\Ui + 02^2 + • • • a r u r eC | a>0 and y, a t- = !}• (2) 

e=i 

We restrict our attention to QCQP problems for which the origin of the complex plane does 
not belong to the relative interior of the convex hull of the set {CV,-, [£?*]#, k G /C} for any edge 
(i,j) in Q(P). For any edge (i, j) in Q(P), consider the points CV, and [Ck]ij,k G /C on the 
complex plane. The convex hull of these points is either a single point, a line segment, or a 
convex polytope. The condition states that (a) if the hull is a single point, then that point is not 
the origin, (b) if the hull is a line segment, then the origin is either an extreme point of the line 
segment or the origin does not lie on the line segment, and (c) if the hull is a convex polytope 
then the origin is either outside or on the boundary of this polytope. This is illustrated in Figure 

m 

We also limit the discussion to QCQPs for which the set of feasible solutions is bounded and 
has a strictly feasible point, i.e., there exists ieC™ such that x n CkX < bk for all k G fC. 
To summarize, consider QCQPs that satisfy the following. 
Condition 1: (a) G(P) is connected and acyclic. 

(b) For any edge (i,j) in G(P), the origin is not in the relative interior of the convex hull of 
{Cii,[Cfe]o-,A;e/C}. 

(c) The set of feasible solutions of P is bounded and has a strictly feasible point. 
The main result of the paper is the following theorem and its application. 
Theorem 1: All QCQPs P that satisfy condition [T] can be solved in polynomial time. 

For a continuous optimization problem, we say it can be solved in polynomial time if given 
any ( > 0, there is an algorithm that finds a feasible solution to the optimization problem 
with an objective value within £ of the theoretical optimum in polynomial time J5j, [35 1, [36|. 



For QCQPs P that satisfy condition [T] Theorem [j] says that we can construct such a point in 
polynomial time. 

To solve P, we use its convex relaxation that can be solved in polynomial time. The relaxation 
is said to be exact, if there exists an optimal solution of the relaxation that can be mapped to an 
optimal solution of P. An exact (convex) relaxation by itself does not however guarantee that P 
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(a) Origin does not belong to the relative interior of the convex hull of these points. 




- - *■ 

Re 




(b) Origin lies in the relative interior of the convex hull of these points. 



Fig. 1: Convex hull of {CV,, [C^Jy, k G JC\ in condition 
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can be solved in polynomial time. This is the case when the set of optimizers of the relaxation 
contains solutions that cannot be mapped to a feasible point in P and there may or may not be 
a polynomial time algorithm to find a correct optimum among that set of optimizers [j] Theorem 
[T] asserts that when condition 1 holds, not only the convex relaxation of P is exact, but it can 
also be solved in polynomial time. See Remark 1 at the end of this section for more details. 

In the remainder of this section, we prove Theorem [T] for the case where C is positive definite. 
The proof for the positive semidefinite case is presented in the appendix. Our proof requires the 
following result on Hermitian matrices, that is presented here and proved in the appendix. 

Lemma 2: If H \ >z and H 2 >z are two n x n matrices, 

tx{H X H 2 ) > Pmin[#l] Pmax[# 2 ]. 

where p m i n [H] and p max [H] respectively denote the minimum and maximum eigenvalues of any 
Hermitian matrix H. 

A. Proof of Theorem \l\for C >- 0: 

For C y 0, we prove the result more generally by relaxing the condition that the feasible 
region of P is bounded. Consider the following semidefinite program RP where W is an n x n 
complex positive semidefinite matrix. 
Relaxed Problem RP: 

minimize tr(CW) 
wyo 

subject to: tr (C k W) <b k , ke hC. (3) 

RP is a convex relaxation of P (5J, p6| . Define p* and r* as the optimum values of the objective 
functions for problems P and RP respectively. 

Lemma 3: p*, r* are finite and p* > r*. If W* solves RP optimally and rank W 7 * < 1, then 
p* = r* and problem P has an exact SDR. 

Proof: The objective functions of P and RP are nonnegative and hence p* and r* are 
finite. Given any feasible solution x of P, W := xx n is a feasible solution of RP. Hence RP 

1 If every optimal solution of the relaxation can be mapped to a solution of P, then clearly P can be solved in polynomial 
time. This is sufficient but not a necessary condition for a polynomial time solution. 
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is feasible and p* > r*. If rank W* = 0, then W* = 0, and an optimal solution to P is x* = 0, 
and therefore r* = p*. If rank W 7 * = 1 then W* has a unique decomposition W 7 * = x^x^-, where 

r* = tr(CW*) = a^Ca;* = p*. 

■ 
Next, we show that there exists a finite W* that solves RP optimally and has rank W* < 1. 

Let the Lagrange multipliers for the inequalities in (|3]) be A& > for each k E /C. Then the 
Lagrangian dual of i?P is 
Dual problem DP: 



maximize 

A>0 

fce/c 



— 2_^ ^ k ^ k 

keK, 

subject to C + ^ A fe Cfc ^ 0. 
fce/c 

For convenience, we introduce the n x n matrix A(\) defined as: 

A(\):=C + J2^C k , (4) 

kdK 

Define a function T from Hermitian matrices to undirected graphs as follows. For any n x n 
Hermitian matrix H, the graph J-'(H) on the vertex set [n] satisfies 

(i, j) is an edge in F(H) •<=>- i ^ j and Hij ^ 0. (5) 

From the definitions of J 7 , Q and A(X), it follows that for any A, T(A(X)) is a subgraph of 
£(-P). For some values of A however, edge (i, j) may exist in G(P) but not in ^(^(A)); in this 
case ^(^(A)) is acyclic but may not be connected, and hence it may be a forest of two or more 
disconnected trees rather than a single connected tree that spans all vertices in the graph. Now, 
we present a lemma about the connectedness of the graph ^(^(A)). 

Lemma 4: For all A ^> 0, J r (A(A)) is connected. 

Proof: Consider any edge (i, j) in Q(P). From condition [TJ the origin is not in the relative 
interior of the convex hull of {CV,, [Cj.]^-, k E /C}. Using Q, we have [A(A)]jj ^ and hence 
(2, j) is an edge of ^(^(A)). Thus ^(^(A)) is identical to Q(P) which is a tree that spans all 
the vertices of the graph. ■ 

Next we characterize the relationship between the optimal points of RP and DP. Let d* 
denote the optimal value of the objective function of DP. 
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Lemma 5: r* = d* and RP/DP has a finite primal dual optimal point (W 7 *, A*). 

Proof: To prove this, we first show that DP is strictly feasible. At A = 0, p m in[/4(A)] = 
Pmin[C] > 0. For a sufficiently small A' 3> 0, p m - m [A(\')} > and hence A(X') y 0. This implies 
A' is a strictly feasible point of DP. Also since P is assumed to be strictly feasible, RP is 
strictly feasible and it has a finite optimum r*. The rest follows from Slater's condition p6| . ■ 
For convenience, define A* := /L(A*). Now we turn our attention to the graph of A*, i.e., F^A*) 
to further analyze the primal dual optimum point (W*, A*) of RP/DP. 

Lemma 6: If J 7 (A*) is connected then rank W* < 1. 



Proof: We observe that rank A* > n— 1. This follows from a result in the literature [37|, [38 
Theorem 3.4] and |39} Corollary 3.9] that states that for any n x n positive semidefinite matrix 
H where the associated graph ^(H) is a connected acyclic graph (i.e., a tree), rank H > n— 1. 
Next we show that rank H 7 * < 1. The complementary slackness condition for optimality of 
(W*, A*) implies 

tr(,W*) = 0. 

Let W 7 * = J2i PiWiwf 1 be the spectral decomposition of W*. Then, 

tr(A»W*) = Yl pi w ? A * w i = °- 

j 

Since A* ^ 0, the eigenvectors iUj of W 7 * corresponding to nonzero eigenvalues pi are all in the 
null space of A*. The rank of A* is at least n — 1 and hence its null space has dimension at 
most 1, from which it follows that rank W* < 1. ■ 

J- (A*) can be connected in one of two ways: (a) the origin of the complex plane lies strictly 
outside the convex hull of the set of points {CV,, [Ck]ij ^ 0, k E fC} for all edges (i,j) in G(P), 
or (b) A* ^> (from lemma [4]). In both cases, lemma [6] guarantees that rank W 7 * < 1. 

If the origin lies on the boundary of the convex hull however, then J- (A*) may not be connected 
when A* is not element-wise strictly positive and therefore rank W* < 1 may not hold. In that 
case, we cannot obtain an optimum solution of P from the optimum solution of RP. We deal 



with this case where J- (A*) is disconnected by using a perturbation [40 1, [41 1 of RP/DP so 
that J 7 (A*) is connected in the perturbed problem. Then we recover an optimal solution W* for 
RP such that rank W* < 1 from the perturbed problem. 
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Define the perturbed problems for parameter e > 0: 
Perturbed relaxed problem RP £ : 



minimize tr(CW) - e V* \b k - tr (C k W)\ 
W)-0 ^— ' 

fc6/C 

subject to: tr (C fc W) < 6 fe , fee /C. 



Perturbed dual problem DP £ 



maximize 

A 



-^Xkbk 



keK. 



subject to A (A) y 0, \ k > e, k € tC. 

For any variable z in the original problem, let z £ denote the corresponding variable in the 
perturbed problem with perturbation parameter e. 

Lemma 7: There exists a e > 0, such that for all e in (0,e ): 

1) r £ = d% and RP £ / DP £ has a finite primal dual optimal point (W £ , A^). 

2) J-'(^) is connected and rank W £ < 1. 

Proof: We choose £ as follows. For e > 0, observe that 

fce/c 
where 1 is a vector of all ones of appropriate size. Then p min [A(£l)] is a continuous function 
of e that has the value p m m[C] > at e = 0. Choose £ sufficiently small such that p min [A(el)] 
is strictly bounded away from 0, i.e., 

min p m[n [A(el)] > 0. 

ee[0,eo] 

Consider any e in (0, Eq). The feasible sets of RP and RP £ are identical. Since A(el) y 
and W y for a feasible point W of RP (and RP £ ), 



fce/c 



tr(CW) - e ^ [6 fc - tr (C fc W^)] = tr [A(el) W] - £^ 6 fc > -e 

fcex: fceK 

and hence r^ is finite. RP £ / DP £ are strictly feasible and r^ is bounded below. The first part of 
lemma [7] then follows from Slater's condition (36). 



To prove the second part of lemma [7J note that A* > el ^> 0. Lemmas [4] and [6] applied to 
RP £ proves the claim. 
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We have shown that for all e in a nonempty interval (0, e ), the optimal solution Wl of RP £ 
has rank at most 1. Now we analyze the behavior of Wl as e converges to zero. Let {e(}^L 1 be a 
decreasing sequence such that Ei — > as £ — > oo. Consider the sequence of matrices {Wl 1 }^; 
every matrix in this sequence has rank at most 1. In the next lemma we show that this sequence 
has a convergent subsequence and the limit of this subsequence solves RP optimally. 

Lemma 8: {Wl 1 } ( ^ =1 has a convergent subsequence. The limit point W of this subsequence 
solves RP optimally and satisfies rank W < 1. 

Proof: Consider any e in (0, E ). We first show that Wl is bounded, independent of e. For 
any W in the feasible set of RP (and RP 6 ), 



fce/c 



tr(CW) - e V [b k - tr (C fc W)] < tr(CW). 



>o 



Minimizing both sides over the feasible set of i?P (and RP £ ), we have 

rf < r*, 



(6) 



that implies 



ti[A(el)W; 



+ E^b k 



keK 



E»» 



fce/c 



< r* + £ 
Also, Lemma [2] implies 

tr[A(el) W:\ > Pmin [A(El)} p max [W* 
min p min [A(E'l)} 

e'e[0,£ ] 



(7) 



Pmax[^]. 



(8) 



>0 by construction. 



Combining equations ^ and ([8]), we obtain a bound for pmax[WJ]» independent of e: 

^maxL"* J ^ . r /!/ /-i \1 ' 

min^o^/w^e'l)] 
Thus WJ is bounded and rank Wl < 1. Since the set of positive semidefinite matrices with 



rank at most 1 is closed [42 1, Wl lies in a compact set and {W^}^ has a convergent subse- 
quence. Let the limit of this subsequence be W . Then W is feasible for RP and rank W < 1. 
Next, we prove that r* = tr(CW) and hence W solves RP optimally. 
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From (|6]), 

tr(CWf ) - e J2 [h ~ tr {C k W*)\ < r„. (9) 






fce/c 
Taking limit over the convergent subsequence of {W^}^, we get tr(CW) < r*. Also, r* is 
the optimum value of PP and hence r* < tr(CW). This completes the proof of lemma |8| ■ 

So far we have shown that RP has a minimizer W that satisfies rank W < 1 and p* = r*, 
i.e., SDR of P is exact. But it is, in general, hard to guarantee that solving RP would yield the 
minimum rank optimizer if the set of optimizers of RP is non-unique. In that case, RP cannot 
be directly used to solve P in polynomial time. In what follows, we present an algorithm to 
solve P in polynomial time. 

First, solve RP in polynomial time to obtain r*. If the associated optimizer W* has rank at 
most 1, then construct s* from W* as in lemma [3] We have then found x* in polynomial time 
that solves P optimally. If however rank W* > 1, choose e as given above and solve RP e ° in 
polynomial time. For any e in (0,£ ), 

rl = tr(CWf ) - e ^ [b k - tr(C fc W^)] < r„ < tr(C^ e ), (10) 

keK. 

where the first inequality follows from ([6]) and the second one follows from the fact that r* is 
the optimum value of RP. Also, comparing the objective function values of RP £ and RP £ ° at 
Wl and W £ ° respectively, we have 

J2 [h - tr(C k W:)} <J2^- tr(C k W:°)] . (11) 

keK. keK. 



Combining ( |T0| ) and ( [TT] ), we have 

|r # - tr(C^ £ )| < sJ2 lh ~ tr(CWf)] • 

keK. 

Given any ( > 0, choose e in (0,£:o) such that Y.k&?c [h ~ tr(C fc WJ )] < (. Now solve i?P e in 
polynomial time to get W £ that satisfies rank W £ < 1 and compute x% from it. Then x% is a 
feasible point of P and 

P* < (xl) H C(xl) < P * + (. 

Also, we have computed x% in polynomial time. This completes the proof of Theorem [T] for the 
case where C is positive definite. We remark that a perturbation by an arbitrary small e > 
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can be represented by treating each perturbed scalar variable as a pair [a, a'} to represent a + ea' 
when solving RP using any standard polynomial time algorithm like the interior-point method 
(TJ-[|3j, where the pairs [a, a'} are ordered lexicographically [43 1. 

The proof extending the theorem to the case where C is positive semidefinite is given in the 
appendix. 

Remark 1: The strict feasibility of P in Condition [T] is required to solve P in polynomial 
time. If we relax that constraint, it can be shown that there still exists a positive semidefinite 
matrix W* that solves RP optimally and rank W* < 1, i.e., P has an exact SDR. There might 
also be other optimal solutions of RP that do not satisfy the rank condition and hence cannot 
be mapped to an optimal solution of P. Solving for a low-rank optimizer arbitrarily closely in 
polynomial time is hard to guarantee and is a direction for future work. 



III. Optimal Power Flow: An application 

In this section, we apply the results of Section [n] to the optimal power flow (OPF) problem. 
We start by summarizing some of the recent results on OPF relaxations in Section III-A| In 



Section |III-B| we formulate OPF as a QCQP In Section |III-C| we restrict our attention to OPF 
over radial networks, which are the networks commonly found in distribution circuits, and use 
Theorem [T] to characterize a set of conditions under which OPF can be solved efficiently. 



A. Prior work 

As previously discussed, OPF can be cast as a QCQP. Various non-linear programming 



techniques have been applied to the resulting nonconvex problem, e.g., in [44 1— [46]. An SDR 



for OPF has been explored in [29|, [30| and their simulations indicate that the SDR provides 



an exact solution of original OPF for many of the IEEE test systems [ |47[ |. The authors in pT| , 



1 48 1 propose to solve the convex Lagrangian dual of OPF and derive a sufficient condition under 
which an optimal solution of OPF can be recovered from an optimal dual solution. Though SDR 
recovers a solution to OPF on IEEE test systems, it does not work on all problem instances. 



Such limitations have been most recently reported in [49|, though the nonconvexity of power 



flow solutions have been studied much earlier, e.g., in p3| , [|50|-[p2|. 

Recently, a series of work has explored a class of problems where such limitations do not 



apply due to the network topology. It has been independently reported in [32|-[34| that the SDR 
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of OPF is exact for radial networks provided certain conditions on the power flow constraints 
are satisfied. A different approach to OPF has been explored using the branch flow model, 



first introduced in [53 1, [54|. While [55| studies a linear approximation of this model, various 



relaxations based on second-order cone programming (SOCP) have been proposed in [56|-[59|. 



Authors in [57|-[59| prove that this relaxation is exact for radial networks when there are no 
upper bounds on loads, or when there are no upper bounds on voltage magnitudes. 



Motivated by the results in [57], a more general branch flow model is introduced in [ |60[ 
for the power flow analysis and optimization for both radial and meshed networks. The precise 
relationships between the SOCP relaxations and the SDR for the OPF problem has been recently 
identified in [loTl. 



B. Problem Formulation 

Consider a power system network with n nodes (buses). The admittance-to-ground at bus i 
is yu and the admittance of the line between connected nodes i and j (denoted by i ~ j) is 
Dij = gij — ibij. We assume both g^ > and b t j > 0, i.e., the lines are resistive and inductive. 
Define the corresponding n x n admittance matrix Y as 



Y; 



Vii + ^Vij, if 



J- 



j~* 



>j 



'Viji 







if i ^ j and i 
otherwise. 



(12) 



Remark 2: Y is symmetric but not necessarily Hermitian. 

The remaining circuit parameters and their relations are defined as follows. 

• V and I are n-dimensional complex voltage and current vectors, where V&, h denote the 
voltage and the injection current at bus k E [n] respectively. The voltage magnitude at each 
bus is bounded as 

0<W k <\V k \ 2 <W k , ke[n}. 

• S = P + iQ is the n-dimensional complex power vector, where P and Q respectively denote 
the real and reactive powers and 

S k = P k + \Q k = V k I k H , ke[n}. (13) 
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• P k and Q k axe the real and reactive power demands at bus k G [n], which are assumed to 
be fixed and given. 

• P k and Q k are the real and reactive power generation at bus k. They are decision variables 
that satisfy the constraints Pjf < Pg < P k and Q° k < Q k <Q k . 

Power balance at each bus k e [n] requires P k = P k + P k and Q k = Q k + Q k , which leads 
us to define 



P G — P 



D 



=;G 



p k := p; - P k D 

Q_ k := Q_ k ~ Qk 5 Qk := Qk ~~ Qk ■ 



^i 



\D 



The power injections at each bus k E [n] are then bounded as 

P k <Pk<P k , Q k <Qk<Q k - 

The branch power flows and their limits are defined as follows. 



• bj. 



Pij + iQij is the sending-end complex power flow from node i to node j, where 



Pij and Qij are the real and reactive power flows respectively. The real power flows are 
constrained as \Pij\ < Fij where Fy is the line-flow limit between nodes i and j and 



. L 



i.i 



F ■ 

Pij + Pji is the power loss over the line between nodes i and j, satisfying Ly < Li 



where Ly is the thermal line limit and Lij 



La;. Also, observe that since L, 7 > 0, we 



j ji 



have | Pj | < F^, \Pji\ < F j{ if and only if Py < Fij,Pji < F^. 
Let J k = e^e^ where e k is the £>th canonical basis vector in C n . Define Y k 
Substituting these expressions into ( [T3] ) yields 



e k eZY. 



e%VI n e k = tr (VV H (Y H e k e%)) = V n Y k H V 
( \ ( 



Y 



H 



Y? + Y k 



Y 



\ 



+ i 



/ 



V 



H 



y k H - v k 

2i 



V 



(14) 



/ 



where $ fc and ^ k are Hermitian matrices. Thus, the two quantities V n Q k V and V n ^ k V are 
real numbers; moreover 

P k = V H $ k V, Q k = V H V k V. 



January 1, 2013 



DRAFT 



15 



The real power flow from i to j can be expressed as a quadratic form as follows. 

Pij = ReiViiVi - Vj) n y%} = V n M«V 7 (15) 

where M !J is an n x n Hermitian matrix. Further details of the OPF problem formulation are 
provided in the appendix. 

The thermal loss of the line connecting buses i and j is 

L tj = Lji = P l3 + Pji = V H T ij V (16) 

where T^ = T>* : = M** + M* h 0. 

For a Hermitian positive semidefinite n x n matrix C, we have 
Optimal power flow problem OPF: 

minimize V H CV 

yec™ 

subject to: 



P k <V H $ k V<P k , ke[n], (17a) 



Q h < V H ^ k V <Q k} ke [n], (17b) 



-k 



W k <V n J k V<W k , ke[n], (17c) 

V n M ij V < Ftj, i ~ j, (17d) 

V n T ij V<L ij} i~j, (17e) 



where (17a)-(17e) are respectively constraints on the real and reactive powers, the voltage 



magnitudes, the line flows and thermal losses. Note that since T lJ >z 0, ( [16] ) implies that 
Pij + Pji > 0. This means that ( |17d| ) holds if and only if |Py| < Fy, i.e., ( |17d| ) bounds 
the line flows on both ends. 

We do not include line-flow constraints that impose an upper bound on the apparent power 
\ Pfj + Qij on each branch i ~ j. These constraints are not quadratic in voltages and hence 
beyond the scope of our model. 

Remark 3 (Objective Functions): We consider different optimality criteria by changing C as 
follows: 

• Voltages: C = l nxn (identity matrix) where we aim to minimize || V]| 2 = Xlfcl^fel 2 - 
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• Power loss: C = (Y + Y H )/2 where we aim to minimize £\ gu\ V^| 2 + J2i<j 9ij\Vi ~ Vj\ 2 - 

• Production costs: C = J2k c ^k where we aim to minimize Yl<k c kPk' c k ^ 0- 
We assume C is positive semidefinite. 

C. Semidefinite relaxation of OPF over radial networks 

Assume hereafter that OPF is feasible. To conform to the notations of Section [III replace the 



constraint in ( |17a[ ) by the equivalent constraints 

V H [$ k ]V < P k , ke[n] 
V n [-$ k ]V < -P k , ke[n] 



Similarly rewrite ( |17b| ) and ( 17c| ). Then the set of matrices {C k , k G K.} and the set of scalars 
{b k , k G /C} in OPF are defined as 

{C k ,ke1C}:={Q k ,-& k ,* k ,-* k ,J k ,-J k , kE[n}} \J {M ij ,T ij , i ~ j} (18) 

{h,kelC}:={p k ,-P k ,Q k ,-Q k ,W k ,-W k , ke[n}] \J {T j ,V j , i~j} (19) 

We limit the discussion to OPF instances where the graph of the power network is a tree T 
on n nodes. It can be checked that the graph of the problem OPF satisfies 

g{OPF) = r. (20) 

Thus condition [Tfla]) holds for OPF over T. In general, condition [TJ|b]) does not hold for OPF 
over T . To illustrate this point, consider an edge (i,j) in T. The admittance of the line joining 
buses i and j is g^ — ibij. Then [C^]^, k e K, are given as (details are in the appendix): 



(a) 


Ma = 


-9ij/2 + 


ibij/2, 


(b) 


[®j\ij = 


-9ij/ 2 ~ 


ibij/ 2, 


(c) 


[*dy = 


-hi/2 - 


Wj/2, 


(d) 


[*j]« = 


-b tj /2 + 


Wij/2, 


(e) 


[M*fo = 


= -gij/2 + ibij/2, 


(f) 


[M*fo = 


~- -9ij/2 - 


- ibij/2, 


(g) 


[T ij h = 


[T ji h = 


~9ij- 


For 


the objective functions considered, we have 


(a) 


Voltages 


Oij — (J, 
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T J ,T j \C 



^3 



$i, \M 13 



-6/2 



-5/2 



$j, M J * 



#,- 



5 /2... 

6/2 



-*< 



6/2 



-^ 



-6/2- 
~5/2 



5/2 



$z 



^. 



#e 



Fig. 2: Cy and non-zero ([Cfc]^, k E /C) on the complex plane for OPF for a fixed line (i, j} 
in tree T. 



(b) Power loss: Cy = —gij, 

(c) Production costs: C y - = — 5jj(q + Cj)/2 + i6 w -(cj — Cj)/2. 

For the purpose of this discussion, consider power-loss minimization as the objective, i.e., C y - = 
— <7ij. Also, assume ^ > 6 y - > 0. We plot the non-zero (i,j)-th entries of the matrices C^, k E /C 
and C, on the complex plane in Figure [2] and label each point with its corresponding matrix. 
Clearly if we consider all the points [(?&]#, k E /C and Cy, the origin of the complex plane lies 
in the relative interior of the convex hull of these points, i.e., condition [TJb]) does not hold. 

To apply Theorem [l] to OPF, consider the index-set K. C /C so that condition |ljbj) holds for 
Ck,k E K, and C. This corresponds to removing certain inequalities in OPF, i.e., bk = +oo 
for k E K \ K. For example, removing — $, from the set {Ck,k E /C} corresponds to setting 



P; 



-oo. 



Condition [TJjc]) requires that the feasible set of OPF be bounded. This is always the case 



January 1, 2013 



DRAFT 



18 



when W k is finite for all buses fce [n\. We discuss the technical condition of strict feasibility 
in remark HI 

Theorem 9: For IC C /C, suppose condition [l] holds for OPF with C k ,k E K. and C. Then 
OPF can be solved in polynomial time. 

Remark 4: When OPF is strictly feasible, it can be solved in polynomial time. When strict 
feasibility does not hold, it follows from remark [T] that OPF still has an exact SDR but we do 
not guarantee a polynomial time solution. 

We explore, through examples, some constraint patterns for OPF over radial networks where 
a polynomial time solution (for strictly feasible OPF instances) or an exact SDR (for OPFs that 
may not be strictly feasible) is guaranteed. 

Example 1: In Figure [2j consider the (i, j)-th elements of the following set of matrices: 

{$„ $,, %, * j: -%, M { i, M ji , Ti = T*\ C) . 

The origin of the complex plane lies on the boundary of the convex hull of these points but not 
in its relative interior. With this set of points, associate a constraint pattern defined as follows. 
For any point in the diagram that is not a part of this set, the inequality associated with that 
matrix is removed from OPF. For example, the matrices — $-,-, — $» and — \Pj do not feature on 
the list of points. Hence, 

Pj =Ei = % = -oo. (21) 

This can be generalized to a constraint pattern over T by removing the lower bounds on real 
powers at all nodes and the lower bounds on reactive powers at alternate nodes. 

Example 2: Suppose P_ k = Q = — oo for all nodes k in T. This corresponds to considering 
points only on the left-half plane in Figure [2] for all edges (i,j) in T. Clearly, condition [T|b]) 
holds in this case. In Figure [2j we assume gij > bij > 0. Regardless of the ordering between g^ 
and b^ for edges (i, j) in T, the set of points considered in this constraint pattern always lies 
in the left half of the complex plane. 

Removing the lower bounds in real and reactive power can be interpreted as load over- 
satisfaction, i.e., the real and reactive powers supplied to a node k can be greater than their 
respective real and reactive power demands P® and Q®. OPF on a radial network with load 



over-satisfaction can be solved efficiently. This result has been reported in [32]— [34| 
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Example 3: Consider voltage minimization, i.e., C = X nxn . In Figure [2j consider the (i,j)-th 
entries of the following set of matrices: 

1 =^t ) ^ j i ^ j i ^ i i i ' ^ / • 

The origin of the complex plane is again on the boundary of the convex hull of these points. 
The constraint pattern associated with this set of points is 



Pi = Qj = Lij = Lji = Fy = Fji = +oo, and Q. = -oo. 

This constraint pattern is consistent with condition [THb]) over the edge (i,j) and we can construct 
a constraint pattern for the OPF problem. 

IV. Numerical examples 

A. Numerical techniques 

In Section |n} we have identified conditions under which an SDR of a QCQP over an acyclic 
connected graph can be used to solve the nonconvex problem P efficiently. When these conditions 
are not satisfied, the QCQP may not be polynomial time computable or have an exact SDR, i.e. 
RP yields an optimal W* with rank W 7 * > 1. In that case, W 7 * cannot be mapped to an optimal 
a;* for the problem P. In this section we propose a method to construct a feasible solution x 
for P using such an optimal W* of RP. The following relation characterizes the relationship 
between the optimal solution of P and its value at x: 

objective value of RP at W* < optimum objective value of P 

< objective value of P at x. (22) 

In many practical problems where rank H 7 * > 1, the principal eigenvalue of W* is orders of 
magnitude greater than the other eigenvalues. We use the principal eigenvector to search for a 
"nearby" feasible point of P as follows. Let w* G C n be the principal eigenvector of W* and 
define the starting point x of the algorithm as 

x := w*\/tT(CW*). 

This scaling ensures that the objective value at x is the optimum objective value of RP at W*. 
If xq satisfies all constraints in P, then the algorithm ends with x = xq. Otherwise, we construct 
a sequence of points (xi,X2, . . .) where x m+ i is constructed from x m as follows. 

January 1, 2013 DRAFT 



20 



1) For k 6 tC, linearize the function fk(x) = x n C k x around the point x m and call this function 



fi m \x), i.e. 



2) For k E JC, define 



/* ( X ) = X m C kX m + 2 Re [x™C k (x - Xr. 



'h-fi m) (x), if fJT\x)<b_ k , 



4 m) W := 



ifb k <fi m \x)<h, 



fJT\x)-b k , ifb k <fJT\x) 



We can interpret s k (x) as the amount by which the linearized function f k violates the 



inequality constraint b k < f k (x) < bi 



3) Compute x m+1 using 

x m+ i = argmin ^2[s ( k m) \x)} 2 

xe£n k&K, 

subject to: ||x — a; m ||i < 7, 

where ||.||i denotes the i 1 norm and 7 is the maximum allowable step-size. This is a 
parameter for the algorithm and should be chosen such that the linearization f r k n (x) is a 
reasonably good approximation of the quadratic form f k (x) for all k E /C in the i 1 ball 
centered around x m with radius 7. 

4) If x m+ \ satisfies all constraints in P, then the algorithm ends with x = x m+ i. 

This heuristic approach either ends at a feasible point x of P or it fails to produce one within a 
fixed number of iterations. In the next section, we show that this technique performs quite well 
for numerical OPF examples where the SDR yields an optimal solution W* with rank more 
than 1. 

B. OPF test examples 



The SDR of OPF and the techniques described in section III are illustrated on a sample 



distribution circuit from Southern California and randomly generated radial circuits. The semidef- 



inite program is solved in MATLAB using YALMIP [62|. If the solution yielded W* such that 
rank H 7 * < 1 then the optimal voltage profile (K) to the OPF problem is calculated from 
W* = (K)(K) W - If the optimal W 7 * does not satisfy the rank condition, the heuristic approach 
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described above is used to find a feasible point of OPF. The feasible point obtained may not be 
optimal for the original problem, so we characterized its sub-optimality by defining the following 

quantity. 

Objective value at heuristically reached feasible point 

77 := 1. 

Objective value at optimal point of relaxed problem 

Smaller values of 77 indicate that the feasible point obtained using the algorithm is close to the 
theoretical optimum of OPF. 

Throughout this section, let y = (a, b) denote a y drawn from a uniform distribution over the 
interval [a,b]. Using this notation, we describe the test systems used for simulations. 

1) SoCal Distribution Circuit: The sample industrial distribution system in Southern California 



has been previously reported in [57|. It has a peak load of approximately 11.3 MW and 
installed PV generation capacity of 6.4MW. We modified this circuit by removing the 
30MW load at the substation bus (that represented other distribution circuits fed by the 
same substation) and simulated it with the parameters provided in Table [TJ To scale the 
problem correctly, the problem was cast in per unit (p.u.) quantities using base values given 
in Table E 
2) Random Test Circuits: These circuits are generated using parameters typical of sparsely 
loaded rural circuits, as detailed in [63] and employed (with suitable modifications) in [64], 



[65]. Around 15-60% of the nodes are assumed to have 2 kW of PV capacity. The remaining 

parameters of these systems are described in Table [TJ 
The tests are run with both voltage and power-loss minimization as objective functions. The 
optimization results are summarized in Table [TTJ. For power-loss minimization, we always obtain 
a rank 1 optimal W*. 

For voltage minimization, however, we obtain optimal solutions that violate the rank condition. 
In these cases, the heuristic approach is used to find a feasible point of OPF. We construct the 
solution based on the complex voltage V k = \Vk\e ldk at bus k E [n]. For the heuristic algorithm, 
define 

x:=(\v 2 \,\v 3 \,...,\v n \,e 2 ,e 3 ,...,e n ), 

and set the parameter 7 = +00. In the examples studied, this approach always yields a feasible 
point within 5 iterations. From Table [IlJ the values obtained for r/ indicate that our algorithm 
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Test system 


SoCal distribution circuit 


Random radial networks 


Number of nodes (n) 


47 


50-150 


Line impedances (y%j) 


H 


Table 1] 


(0.33 + i0.38)fi/Jfem, length = (0.2km, 0.3km) 


Voltage limits \/Wk, \/W^ 


1 ± 0.05 p.u. at all nodes. 


1 ± 0.05 p.u. at all nodes. 


Real power demand P® 




57 


Table 1] 


(0,4.5fcW) 


Reactive power demand Qj? 


Computed with p.f. = (0.80,0.98) lagging 


(0.2P, D ,0.3P, D ) 


Real power gen. limits 

~P G P G 
Pk > ilk 


q 

PV nodes: P k = (0.2, 1.0) times capacity, 

Substation node: P° = WMW. 

At all nodes, P° = 0. 


PV nodes: P^ = (0,2feW), 

/-* 

Substation node: P 1 scaled with n. 

At all nodes, P% = 0. 


Reactive power gen. limits 


Q G = 0.3P^, Q^ = -0.3P^ at all nodes. 


Q°k = 0.3P°, Q° = -0.3P^ at all nodes. 


Base quantities 


Phase = 1MW, Vbase = 12.35fcV(L - L) . 


P base = 1MW, V b ase = 12A7kV(L - L). 



TABLE I: Circuit Parameters for SDR of OPF 



Test system 


SoCal distribution circuit 


Random radial networks 


Minimize 


Power-loss 


Voltage 


Power-loss 


Voltage 


rank W* 


1 


> 1 


1 


> 1 


Mean 77 


N/A 


1.8% 


N/A 


0.5% 


Maximum r) 


N/A 


4.1% 


N/A 


1.5% 



TABLE II: Summary of simulation results 

finds a feasible point of OPF with an objective value close to the theoretical optimum and 
hence performs well. A general bound on the performance of this heuristic technique, however, 
remains an open question. 

V. Conclusion 

QCQP problems are generally non-convex and NP-hard. This paper proves that a certain class 
of QCQP problems are solvable in polynomial-time. We have applied this result to the optimal 
power flow problem and derived a set of conditions under which this nonconvex problem admits 
an efficient solution. For problems that do not satisfy our sufficient conditions, we provide a 
heuristic technique to find a feasible solution. Simulations suggest that this method often finds 
a near-optimal solution for the OPF problem. 
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VI. Appendix 
A. Proof of lemma [2]" 

The result is restated here for convenience: If H\ >z and H 2 b are two n x n matrices, 

tr(ififf 2 ) > p min [^l] Pmax[^2]- 

Proof: Let the spectral decomposition of i?2 be 



# 2 = $> 

»=i 

where ||uj|| = 1 for all 1 < i < n. Then: 



1-L 



tr(H 1 H 2 ) = J2pi {ufH x Ui) 

8=1 

n 

> ^Pi Pmin[#i] 
i=l 

> Pmin^l] Pmaxf^]- 



5. Proof of Theorem \l\for C >z 0: 

The sketch of the proof is as follows. Perturb RP so that the matrix in the objective function 
is positive definite. From our previous analysis this perturbed problem has a finite optimizer that 
has rank at most 1. Also, this perturbed problem can be solved in polynomial time. Using the 
solutions from the perturbed problems, we construct an optimum solution of P in polynomial 
time. 

In particular, consider the perturbed problems for 5 > 0: 
Perturbed primal problem P^: 

c H ( 
Mi 



minimize x H (C + 51) x 

xec n 



subject to: x C k x < b k , k E K. 
Perturbed relaxed problem RP^: 

minimize trIYC + 5X)W] 

subject to: tr (C k W) <b k , k G /C, 
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where X is the nxn identity matrix. For any variable z in P/RP, let z^ denote the corresponding 
variable in P^/RP {5) . 

The matrix C + 5X is positive definite for all 5 > 0. There exists W* >z that solves RP^' 
and rank W* < 1. Let the spectral decomposition of Wi be W* = (x* )(#* ) w . From 
Lemma |3j we have 

P? = {x^) H {C + £T)(x?) = rf ) = tr[(C + <SZ)Wf>]. 

The feasible regions of P and P^ are the same and are bounded. Then xl lies in a compact 
space, independent of 5. Taking a sequence 5 — > 0, the corresponding sequence of x\ has a 
convergent subsequence. Let the limit point of this subsequence be x. Then x is feasible for P. 
We now show that x solves P optimally. 

The objective function of P( s > increases with 5. Therefore, 



P* 



<x n Cx<pi 5) =rf ) . (23) 



We wish to show that the first inequality is an equality. Suppose on the contrary p* < x H Cx. 
Let x[ be any finite optimizer of P. Then, 

P* = «) H C(x'J < x n Cx, (24) 

and we can choose a sufficiently small S > 0, such that 

P* + 5«) H «) < x n Cx < rl s \ 

This follows from equations ( |23] ) and ( |24| ). For this 5, [(XXx'J^ + (xi )(xi ) n ]/2 is a feasible 
point of RP^ and satisfies 

r{ s) < tr 



(r + ()T) . «)«) W + (^ ) )(^ ) ) W 



= ^w + i [p. + *oo V. 

This is a contradiction and hence p* = x n Cx. 

Now, we show how to use this perturbation technique to solve P in polynomial time. Solve 
RP to get p* = r*. If the optimizer W* of RP has rank at most 1, compute x* from W* as in 
lemma [3] then we have solved P in polynomial time. If it does not satisfy the rank condition, 
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then choose an arbitrary 5 > and solve RP^ S °> in polynomial time to get the minimizer W* 
and the minimum r* ° —pi '. For any 5 in (0, S ), 



p, = n < tr (CWi S) ) < p 



(S) 



,(*) 



Also, p* is convex in 5 [41 1 and hence 



p^<p* + t (p^ o) -p*)- 



(25) 



(26) 



Given ( > 0, choose 5 sufficiently small so that 4- (pi — p*) < C F° r this 5, solve RP^ 
arbitrarily closely in polynomial time to get W 7 * that has rank at most 1 and compute xi . 



From equations ([25]) and ( |26| ), x* satisfies 

K4 5) m*i 5 ))-j9,|<c. 

Thus a;* is a feasible point of P that has a value of the objective function within £ of the 
theoretical optimum and we have shown a polynomial time algorithm to compute it. 
This completes the proof of Theorem [TJ 



C. Matrices involved in OPF: 



Here we compute the (i,j)-th entries of {Ck,k G /C} for OPF. From ( fT4| ), ( [15] ), ( fT6[ ), we 
have the following relations for k e [n], and (p,q) and (z, j) in T: 



[$ 



fcjij 



[* 



fcjij 



[M pq ] 



1 

2 


Y - 


= §(-#*+*&«) 


if fc = i 


lyH 

2 1 ij 


= |(-Py-i6ii) 


if k = j 







otherwise, 


2i 2 *J 


= l(-bij-i9ij) 


if fc = i 


1 yH 

2i J ij 


= K-bij + igij) 


if k = j 







otherwise 




9pq 


if %-. 


= 3 =P 




!(- 


-g pq + ib pq ) if (i 


J) = (P, ?) 



'./ 



§(-&*- i&w) if (i,j) = (q,p) 
otherwise, 



(27) 



(28) 



(29) 
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9 pq if i = J =P or i=j = q 

l TP % = \-g„ X(i,j) = (p,q)<x(iJ) = (q,p) OO) 

otherwise. 
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